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Abstract 

A thermomechanical model of continuous fluid media based on second gradient 
theory is used to study motions in liquid-vapor interfaces. At equilibrium, the model 
is shown to be equivalent to mean-field molecular theories of capillarity. In such 
fluids, conservative motions verify first integrals that lead to Kelvin circulation 
theorems and potential equations. The dynamical surface tension of liquid-vapor 
interfaces is deduced from viscous fluid equations. The result provides and explains 
the Marangoni effect. 
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1 Introduction 



Usually, liquid-vapor interfaces are represented by material surfaces endowed 
with an energy related to the Laplace surface tension. Several studies con- 
ducted in the fields of fluid mechanics and thermodynamics represent the 
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interface as a surface separating two media and having its own energy den- 



critical temperature, detailed theoretical or experimental studies show that 
the capillary layer thickness ranges over a few molecular beams PI8f9"] . For 
fluid conditions defined by temperature and density, molecular models such as 
those used in gas kinetic theory lead to state laws associated with non-convex 
internal energies as the one associated with the van der Waals model for the 
pressure [9|10|lllfl2fl3fl4j . Such models provide a direct verification of the 
Maxwell rule applied to isothermal liquid- vapor transition [T5]; nonetheless, 
they present two disadvantages: 

First, for density values ranging between the vapor and liquid bulk densities, 
the pressure may become negative; however, traction created by simple phys- 
ical experiments can lead to such pressure values [9|12fl6j . 




Second, in the field between vapor and liquid, the internal energy distribution 
cannot be represented by a convex surface associated with the density and 
entropy and seems to be inconsistent with the existence of equilibrium state. 
To overcome these disadvantages, the thermodynamical investigation uses a 
plane domain in place of the portion corresponding to the non-convex inter- 
nal energy. The fluid can no longer be considered as a continuous medium. 
The interface is represented by a material surface with a null thickness. In 
such the only possible representation of the dynamic behavior of the 

interface is a discontinuous surface and its essential structure remains un- 
known. At equilibrium, the disadvantages can be eliminated by appropriately 
modifying the stress tensor of the capillary layer which is then expressed in 
an anisotropic form. As a consequence, the energy of the continuous medium 
must be modified [T0lfTTlfT8] . 

If the fluid interface must be represented as a continuous medium, how can 
the stress tensor of the capillary layer be written in a dynamic expression [I]? 
In previous articles, we proposed a dynamic theory named internal capillarity 
based on the second gradient of deformation of the medium [T9"|2U|I21] . Elabo- 
rated in the sixties for the static case, the theory points out the advantage of 
using a three-dimensional approach of the capillarity in a continuous mechan- 
ical model [TTf2"2"] . The deduced equations of equilibrium provide a satisfying 
representation of isothermal liquid-vapor equilibrium states. The approach is 
not new, and dates back to van der Waals and Korteweg [23121] • Calm and 
Hilliard again used such an approach in 1959 in reference to the free energy 
|10j . and went back to a Landau- Ginsburg model [14] . The representation pro- 
posed in the present study is based on the notion of internal energy which is 
more useful when the temperature is not uniform. The apparent contradic- 
tion between the Korteweg classical stress theory and the Clausius-Duhem 
inequality makes the dynamic study of phase transformation more difficult 
[25], but the use of the capillary fluid model rectifies this apparent anomaly 
[TTfT9"] . In the more general case of internal capillarity, by representing energy 
in terms of second gradient of deformation [2S] and by applying a simple al- 



sity and characteristic behaviors 




. When working far from the 
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gebraic identity, it is possible to draw a relationship between energy equation, 
motion equation, mass conservation equation and entropy [T9]. For a conserva- 
tive fluid, an additional term with a heat flux dimension must be introduced 
into the energy equation. In the case of viscous fluid, previous results pro- 
vide a set of equations which do not modify the Clausius-Duhem inequality 
and make them compatible with the second law of thermodynamics. In the 
non-viscous case, classical fluid motions and motions of fluids endowed with 
internal capillarity reveal a common structure that induces a thermodynamic 
form of equation of motion and potential equations [T7|27|28j . This leads to 
the same classification of motions, generalization of Kelvin theorems, Crocco- 
Vaszonyi equation and first integrals conserved across the interfaces [291130] . 
By representing the internal energy as a function of entropy, density, density 
gradient and by using a single constant A, the resulting equations are thereby 
identical to those obtained with molecular models of mean field theory in the 
isothermal case [9]. 

For a surface area relatively large with respect to the capillary layer thick- 
ness, the surface tension is calculated by integration throughout the interface. 
Lengthwise, the surface tension is not constant and depends on the dynamical 
distribution of density and temperature. This dynamical distribution based 
on equations of motion, calls upon a Navier-Stokes viscosity. For interfaces 
in isothermal equilibrium, the results are classical PflU|ll2|ll3|IT%] . The study 
of motion in the interface without mass transfer needs the surface tension 
gradient and the velocity gradients associated with the dynamic viscosity. 
The ensuing Marangoni effect has been interpreted by using a limit analysis 
wherein the approximate quantities correspond to the physical dimensions of 
the interface. When the mass flow across the interface is not null, a general 
dynamic form of the Laplace equation is given [21]. The method herein pre- 
sented is completely different from classical calculations based on the balance 
equations established for both sides of a discontinuous surface, which takes 
account of density variations only as a difference across the interface [21] • In 
the particular case of isothermal liquid-vapor equilibria, an invariant integral 
of motions compatible with the interface coincides with a generalization of the 
rule advanced by Maxwell, associated with phase transitions [15] . 

The model of a viscous fluid endowed with internal capillarity is therefore 
substantiated by the consequences verified in both equilibrium states and 
non-isothermal motions. The model provides a better understanding of the 
behavior of liquid-vapor interfaces in motion. It at least gives a partial an- 
swer to the question: Is the fluid at the interface rigid or moving ? |32j and 
proposes a theory which takes account of the stress tensor and dynamics in 
the structure of the interface [J]. The resulting behaviors are not the laws 
governing classical fluids since they include an anisotropic stress tensor in the 
momentum equation and an additional heat flux term in the energy equation. 
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2 Equations of motion for a fluid endowed with internal capillarity 



2.1 Case of conservative fluids 



The second gradient theory conceptually more straightforward than the 
Laplace theory, can be used to construct a theory of capillarity. In the present 
text the only addition with respect to compressible fluids is an internal mass 
energy which is a function of density p, entropy s and also grad p. The spe- 
cific internal energy e characterizes both the compressibility and capillarity 
properties of the fluid [33J, 

£ = f(ft,s,P), 

where 

P = (gradp) 2 . 
The equation of motion for a conservative fluid writes 

p a = div a — p grad Q, (1) 

where a denotes the acceleration vector, Q the extraneous force potential and 
er the general stress tensor 

cr = — p I — A (grad p) (grad p) T (2) 

= ~P $ij ~ A p 4 p d , i,j £ {1, 2, 3} 
\ = 2pe' p 

,2 / 



or 

where 
and 



p = p e' p — div(A grad p) (3) 
It should be noted that: 

e = < (4) 

is the Kelvin temperature expressed as a function of p, s and (3. 

Proof: The equation of motion is obtained by using the virtual work principle. 
The virtual displacement, denoted by 5x, has been defined by J. Serrin ([34], 
p. 145). For a fluid endowed with internal capillarity, the virtual work principle 
iv, section 14) is stated as follows: the motion of a fluid is such that: 



5 (^J p (e + Q,) dv^j — J p&§xdv = 0, 



where V is an arbitrary material volume. The variation of entropy ([34], P- 
148) must satisfied the condition 5s = 0; while taking account of 

/ dp \ 85 p dp d5x. 
\ <9x / dx <9x dx ' 
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we obtain 



°P = 2 "a ^T^T gradp 

V ax ax ax / 



and 

8s = e' p 5p + <5s + e'p <5/3, 

as well as eqs. (14.5) and (14.6) in [31]. For virtual displacement <5x null at 
the edge of V and integration by parts, the integrals calculated on the edge of 
V have a null contribution; we obtain: 



/ (pa — diver + pgradfi) Sxdv = 0, 
Jv 



where cr is expressed by Eq. fl2]). Classical methods of the calculus of variations 
lead to Eq. flTJ. It appears that a single term A accounts for the effects of the 
second gradient deformation in the equation of motion. Scalar A depends on 
p, s and /3 as e does. In a study of surface tension based on gas kinetic theory, 
Rocard obtained the same expression as Eq. (j2J) for the stress tensor but with 
A constant. If A is constant, specific internal energy e writes: 

e(j>,s,p) = a(p,s) + — /3. 

The second gradient term A /3/(2p) is simply added to the specific internal en- 
ergy a(p, s) of the classical compressible fluid. The pressure of the compressible 
fluid is denoted P = p 2 a' p and the temperature is denoted 6 = a' s . 

Consequently, 

p = P- A 

For the pressure P, Rocard uses either the van der Waal pressure 

R® 2 

p = "— p - a " 

or other laws of which he provides a comparison. It should be noted that if A is 
constant, this implies that 8 = and there exits a relationship, independent 
of the second gradient terms, between 9,p and s. 



2.2 Case of viscous fluids 



If the fluid is endowed with viscosity, the equation of motion includes not only 
the stress tensor cr, but also the viscous stress tensor er v : 

cr v = v tr D + 2 p D, 
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where D is the deformation tensor, symmetric gradient of the velocity field u 

Of course in second gradient theory, it would be coherent to add terms ac- 
counting for the influence of higher order derivatives of the velocity field to the 
viscous stress tensor o~ v ; this has not been done in the present case. Equation 
(2) is modified by adding the virtual work of the forces of viscosity and we 
obtain: 

pa = div(cr + cr v ) — pgrad Q. 



3 Equation of energy for viscous fluids endowed with internal cap- 
illarity 



Let e be total volumic energy, e = p (^u 2 + e + , q the heat flux vector, r 

P 

the heat supply and h = e H — the specific enthalpy. We define the following 

P 

quantities: 

M = p a — div(cr + o~ v ) + p grad Q, 
B = p + p div u, 

S = p9 s + divg — r — tr(cr v 'D) : 

de r , . , ,. . , . ,. dQ 

E = — + div (e - cr - cr^juj - div(A pgrad p) + divg - r - p— . 
at at 

Theorem: For an internal energy written as in Eq. (1) and for any motion in 
the fluid, the relation 

E -M-u- (^u 2 + h + Q^j B - S = (5) 

is an identity. 

Proof: In the first member, the dissipative terms q,r and o~ v cancel out; this 
also is the case for the extraneous force potential and the inertia terms. After 
having replaced <r,p and 9 by their respective values in Eqs (3-6), it remains 
to prove that the terms from internal energy e also cancel out. These terms 
include the following expressions: 

a) in E : — h div(p e u) — div(cru) — div(A p grad p) 
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= pe + e(p + pdivu) — div(cr)u — tr (^ (T ~q~ ^ ~ div(Apgradp), 

or finally, 

= p(e' s s + e' p p + e'p/3) + e(p + p div u) - div(cr)u + p div u 

+A (grad p) 7 — grad p — div(A pgrad p); 
ox 

b) in M • u : — div(cr)u; 

c) in ^u 2 + /i + fij5: U + ^\(p + p div u); 

d) in S : pe' s s. 

So, the first member of Eq. fl5]) writes: 

A ■ p dV 
ps'pP + — /3 p + A (grad p) T -^— grad p — div(A pgrad p). 



Finally, by taking Eq. (5) into account, we obtain that 

A^gradp) h p div A (grad p) + A (grad p) — grad p — div( A pgrad p) 

is identically null. 

Corollary 1 : For a conservative motion of a fluid endowed with internal 
capillarity, the conservation of entropy along the trajectories is equivalent to 

de ■ dfl 

— + div[(e - tr)u] - div(A p grad p) - p — = 0, (6) 

which is the equation of conservation of energy. This is derived from Eq. ([5]) 
when o~ v = 0, div q = r = and M = 0, B = 0, S = 0. It leads us to add the 
additional term div(A pgrad p) to the equation of energy. The vector A pgrad p 
has the dimension of a heat flux vector and even occurs in the conservative 
case. 

Corollary 2: For a motion of a viscous fluid endowed with internal capillarity, 
the energy equation 

de ■ dVt 

— + div[(e — cr — cr v )u] — div(Apgradp) + divg — r — p—- = 

C/ v \Jv 

is equivalent to the entropy equation 

p6s + div q - r - tr(o\,D) = 0. (7) 
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Equation (JTJ) corresponds to the classical version of the entropy variation ex- 
pressed by the function of dissipation of viscous stresses ip — tr(er v D). 



4 The Planck inequality and the Clausius - Duhem inequality 



For any motion in a viscous fluid endowed with capillarity, it has been assumed 
that 

pOs + divq-r > 0. (8) 

Inequality ([S]) represents the Planck inequality [33]. Let generally write the 
Fourier principle as 

q ■ grad 9^0. 

Then, the Clausius-Duhem inequality can be directly deduced as 

q r 
ps + div---^0. 

When in motion, the fluid is endowed with internal capillarity, the law of 
thermodynamics leads to the existence of a heat flux vector. This heat flux 
vector introduces an additional term into the classical equation of energy, even 
if the fluid is non-viscous. It is possible to extend these results to continuous 
media with an internal energy that contains gradient terms of deformation. 



5 Transformation of the motion equations for fluids endowed with 
internal capillarity 

In this section, we rewrite Eq. (2) in other forms. 
5.1 Case A constant 

Let us note u = Q — A Ap. Then, Eq. (2) can be written 

p a + grad P + p grad u — . (9) 

This is a representation of a perfect fluid where P is the van der Waals pres- 
sure; term u contains all capillarity terms. 
Equations (3) and (5) yield 

try = -P5ij + A U]- p^ k p tk + pp,kk)Sij - Xp^pA . 
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Then, 

a iij = ~ P ,i ■ X l>l ) ., i! 

i.e. 

div cr = — grad P + A p grad( Ap) . 
We can note that for viscous fluids, Eq. writes 

pa + gradP+pgrad(fi - A Ap) - diver,, = 0. (10) 



5.2 Thermodynamic form of the equation of motion 



Commonly - and not only when A is constant - the equation of motion ([TJ can 
be written in the following form 

a = #grad s — grad(/i + O), (11) 

which is the thermodynamic form of the equation of motion. In the case with- 
out capillarity (e = a(p, s), e'p = 0), Eq. ( ITT]) is well-known ([31], p. 171). The 
equation remains valid for fluids endowed with internal capillarity. 

The verification simply consists on writing out the second member of Eq. (ITT]) . 
Let us denote 

A = 9 grad s — grad h i.e. A{ = 9s ^ — h t i . 
According to Eq. 01]), we get: A» = e' s s^ — + — ^ . Then, 

A i = t's S ,i - i £ 's S ,i + £ 'pP,i + £ 'p (P,kP,k) ti - j^P,i - • ( 12 ) 

Using Eq. (jHJ), Eq. ( fT2T) writes (Tijj/p and consequently A = divcr/p, which 
proves Eq. ( ITTi) . 



6 Generalized Kelvin theorems 



Let J be the circulation of the velocity vector along a closed fluid curve C 
convected by the fluid flow 

J = 6 u T dx, 



and 

dJ 



c 



a dx. 

dt Jc 
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Equation flllD implies f a T dx = f 9 ds and consequently yields the following 

Jc Jc 

theorems which are valid for fluids endowed with internal capillarity: 

Theorem 1 : The velocity circulation along a closed isentropic fluid curve is 
constant. 

Corollary 1 : For a homentropic flow [27], the velocity circulation along a 
closed fluid curve is constant. 

Theorem 2: The velocity circulation along a closed isothermal fluid curve is 
constant. 

Corollary 2 : For an isothermal flow, the velocity circulation along a closed 
fluid curve is constant. 



7 Potential equations for conservative fluids and classification of 
motions 



The results of [T7|27f28|l29] can be applied to conservative motions of fluids en- 
dowed with internal capillarity. To any motion of fluids endowed with internal 
capillarity, there correspond scalar potentials (p, ip, r and x verifying 





1 




~ 2 


T 


= 0, 




= 9, 


X 


= 0, 


s 


= 0. 



u 2 - h - a 



The fluid velocity is given by 

u = grad <f + ip grad s + t grad x, 

together with 



dp 
dt 



+ div(pu) = 0. 



(13) 



(14) 



The equations can be used to classify motions in the same way as in the case 
of non-capillary perfect fluids [27,30]. 
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7.1 Homentropic motions 

This case corresponds to s constant in all the fluid. Eq. (|T4|) writes 

u = grad ip + r grad x, 

which leads to the Cauchy theorem [3lj 

d f rot u\ du rot u 
dt \ p J <9x p 

7.£ Oligotropic motions 
Equation ffl^j) writes 



u = grad (p + tp grad s. 
The flows verify the relation rot u . grad s = 0. Surfaces with equal entropy are 
eddy surfaces. The velocity circulation along a closed isentropic fluid curve is 
null. 

Remark: the form of Eqs. (J9j) and (TTTT) allows to generalize the results of [30J . 



8 Generalized Crocco -Vazsonyi equation 



The equation of energy ()6]) can be also written 

de ( dp \ dfl 

— + div[(e + p)u] - div I A — gradp j - p— = 0, (15) 

with e + p = p^u 2 + H0j and denoting H = ~ u 2 + h + Q , Eq. fTT5]) 
writes 

?§■ + div „Hu = ,f + | + div (a I grad . (16) 

The first member is equal to pH+H(p+p div u). Because of mass conservation, 
the definitive form of the equation of energy writes 

dfl dp f. dp , \ 



By taking account of the identity 

9u 



at 



+ rot A u + grad (- u 2 ^j 
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Equation (ITT]) can be written 



— + rot u A u = 6 grad s — grad H. 



(17) 



If the motion is a steady flow, Eq. (JTTJ) makes it possible to conclude that 

H is constant along the stream lines (H = 0) and Eq. f)16p shows that the 
divergence of the partial energy flux pHu is null. Equation f fTTl) leads to the 
Crocco-Vazsonyi equation generalized to fluids endowed with internal capil- 
larity; in stationary motion, 



where s and H are constant along every stream line. 

Remark: Equations (I13p show that along each trajectory, and consequently 
across interfaces, r, x an d the entropy s are constant scalars. They represent 
first integrals of the motion. Other integrals can obviously be found as the 
Kelvin integrals 



are constant along any fluid curve T. Thanks to the Noether theorem, any law 
of conservation can be represented with an invariance group |36j. The law of 
conservation expressed with the Kelvin theorems corresponding to isentropic 
fluid curves is related to the group of permutations associated with particles 
of equal entropy; the group keeps the equations of motion invariant for both 
classical perfect fluid and fluid endowed with internal capillarity, but also for 
any perfect fluid endowed with an internal energy depending on p, s and their 
gradients of any order. It is even tempting to define a general perfect fluid 
by identifying it with an invariant group or with a continuous medium whose 
motions verify the Kelvin theorems [57] . 



9 Dynamic surface tension of liquid-vapor interfaces 

Far from the critical conditions, the thickness of a liquid-vapor interface is 
very small [9fT3] . Outside the capillary layer, density and its spatial derivatives 
have smooth variations. The density in each phase is reached at points located 
within the immediate vicinity of the layer. In our study, bubbles and drops 
of size of a few molecular beams have not been considered. Surfaces of equal 
density materializing the interface are parallel surfaces and are used to define 
a system of orthogonal coordinates [38]. Notations are those used in [39]. The 
subscript 3 refers to the direction normal to surfaces of equal density. Let us 
denote by e 3 the unit vector of this index oriented in the increasing density 



rot u A u = 6 grad s — grad H, 
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direction. The mass flow through the interface is assumed to be null such that 
surfaces of equal density are material surfaces. In the capillary layer, 



divu = 0. (18) 

Extraneous forces being neglected, Eq. (fl8"j) represents velocities compatible 
with the interface. Equation f lTUj) yields: 

1 OP 1 dAp , , _ 

p a tg + grad 4g P = A p grad t9 Ap + 2p (div D) tg , 

where the subscript tg denotes the tangent component to the surfaces of equal 
density. The different linear sizes of the interface should be taken in consider- 
ation. The capillary layer is measured in Angstrom and the surface curvature 
radii are of non-molecular dimensions. The deduced relations result from a 
limit analysis where the parameter related to the thickness of the capillary 
layer tends to zero. 

The subscripts v and / designating the vapor and the liquid bulks, respectively, 
integration of Eq. along the third coordinate line yields, when we assume 
that a 3 is bounded through the interface such that its integral is negligible, 



p - p v = XpAp- A / Ap-^-dx 3 + 2 p (div D) 3 h 3 dx 3 . (20) 

Jx% UX 3 JxVj 



The partial derivatives of velocity with respect to coordinates x\ and 22 are 
assumed to be bounded. Taking account of Eq. ffTSl . the last term in Eq. fl20l 
is negligible. Furthemore, we have the relation: 

R m h 3 dx 3 h 3 \h 3 dx 3 ) z ' 

where R m denotes the mean curvature of the surfaces of equal density in the 
capillary layer, oriented by e 3 [26J. It follows: 



P - P„ = X 



pAp- i(gradp) 2 ) +2 A- Agradp) 2 h 3 dx 3 . (22) 
z J R m Jxi 



Denoting dn = h 3 dx 3 , we finally obtain: 

P - P v = 2 — — where H = X / (grad p) 2 dn, 

Rm Jnv 

which expresses the Laplace equation for motions. We interpret H as the 
dynamical surface tension. 

Remark: Interpretation of the surface tension at equilibrium. 

The eigenvalues of the internal capillary stress tensor are deduced from Eq. 
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©: 

Ai = — p + A (grad p) 2 is the eigenvalue associated with the plane perpendicular 
to grad p. 

A 2 = — p is the eigenvalue associated with the direction grad p. 

In the system of curvilinear coordinates related to the interface, the stress 
tensor writes 

Ai 
Ai 
A 2 

The equilibrium equation for the plane interface is drawn from Eq. (JT]); ne- 
glecting extraneous forces we get: 



A 2 - -P , 

where P denotes the pressure in the liquid and vapor bulks. The force per 
unit length along the edge of the interface (see Fig. 1) is: 



\ih 3 dx 3 



A (grad p) 2 h 3 dx 3 , 



where L denotes the capillary layer thickness. Obviously, P L is negligible. 



vapor 



Capillary Layer 



liquid 



H 



Fig. 1. Interpretation of the surface tension in liquid-vapor interfaces as a layer of 
fluid endowed with internal capillarity. 



Let us denote 



H 



A(gradp) 2 fr 3 dx 3 , 



the force per unit length; then H represents the surface tension of the plane 
interface at equilibrium. 
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10 Practical calculation of the surface tension 



By integrating Eq. ( IT9|) on the third coordinate line, the viscosity term being 
assumed negligible, we obtain 



PP ^dp)- (23) 



. , dP 86 . 1 , .-| i . , OP dp 

Let us notice that — — — is negligible with respect to — — — — . 

ad ox 3 op 0x3 

Taking Eq. ( I2ip into account, a further integration on the third coordinate 
leads to 



^(gradp) 2 = 2A [** (pad p) 2 h 3 dx 3 + p f^p-dp. (24) 

2 R m Jxl J Pv p l 



By denoting x\ the third coordinate of the surface of equal density p = p i} 

A r x 3 

then for X3 G [2C3 , ar^] , the term 2—— / (grad p) 2 h3 dx% is negligible with 

respect to — (gradp) 2 ; that is the case for p x = - (p v + p t ). For p G [p v ,Pi\, 
Eq. (EU writes 

-(gradp) = p / —dp. 

2 J Pv p l 

We obtain the same result for p G [p { , pj\ and finally, 



/ fPi / r u P - P rPi / r u P - Pi \ 

H = ^ U, r 1 du + 1, r L ~ dp du ) ■ (25) 

Relation ( 125]) expresses the surface tension of a liquid-vapor interface for mo- 
tions compatible with the interface. Here P is not only a function of density 
but also temperature and varies according to the location of the point within 
the interface. The viscosity has explicitly disappeared in the expression of H 
(the expression of H does not explicitly take the viscosity into account). In 
the case of plane interface in isothermal equilibrium, the expression given by 
Rocard [9] is found again. The value of the internal capillarity constant can be 
numerically calculated on the basis of experimental values of H and expres- 
sions of P. By injecting the value of A in Eq. (125]) . we can calculate surface 
tension for any dynamic temperature distribution. 



Remark: Eq. ( 123]) yields 



AAp=--^+ \"- 2 dp 

P Pv J Pv P 
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and we deduce 

™ P , Pv Pi 

— dp= 

Pv P 2 Pv Pi 

which represents an invariant corresponding to motions compatible with the 
interface. This invariant leads to the equations (4-11) of [15] demonstrated in 
the specific case of an interface in isothermal equilibrium. The case of plane 
interface leads to the Maxwell equal area rule PITB] . 



11 Marangoni effect in liquid- vapor interfaces 



Let us denote 



r x 3 

Hj (x 3 ) = / . (grad p) 2 h 3 dx 3 , where j G {v, I}. 

Jx\ 

By using calculations as in Section 10, for example for x 3 G [xg, x l 3 ], we obtain 
H v (x 3 ) = V2X ^ J sju J —^dp du^j . 



By injecting Eq. (1221) into Eq. ( 123|) . we obtain: 

A / 2 \ 

P a tg = 2 g rad t s (g rad P) 2 -g rad t s [^- H j( x 3)J-g^d tg P j + 2p(divD) tg . (26) 

Equation fll8p and conditions on partial derivatives of velocity imply: 

f X 3 1 

/ (div D) tg h 3 dx 3 = [De 3 } x J . 

Jxl 3 

Integration of Eq. (|2"oT) yields: 

A f x 3 2 x l 

/ grad t (gradp) h 3 dx 3 + 2[pDe 3 }J = . 



2 Jxi 



Taking account of 



J ' * grad tg (gmd p) 2 h 3 dx 3 = gmd tg |a J \gr&dp) 2 h 3 rfx 3 | , 



A r 
2 



we can conclude 

grad^ + 2[/i J De 3 ]l = 0. 

If the viscosity stress tensor of the vapor bulk is supposed to be negligible, we 
get: 

gmd tg H + 2[pD l e 3 ] = 0. 
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This is the usual Marangoni condition for free boundary problems [6]. In the 
limit case when the viscosity coefficient is null (as for superffuid helium), the 
problem must be proposed in another way: the momentum associated with the 
interface can no longer be neglected and other physical effects must be taken 
into account We must note there exist other phenomenal presentations of 
the Marangoni effect (for example [10]), but to our knowledge, they all consider 
the interface as a discontinuous surface of the fluid medium. The calculations 
performed in the present text do not call upon the use of any linear approx- 
imations, and only take account of the various physical quantities describing 
liquid-vapor interfaces while working far from the critical temperature. The 
use of second gradient theory for representing interfaces has, of course, been 
considered by several authors. In jTTll 1 1 II 12. . their source is found to be in 
free energy which is more directly useable in case of isothermal flows. 



12 Conclusion and further developments 

As we have seen in this revisited version of [33], the second gradient theory, 
conceptually more straightforward than the Laplace theory can be used to 
build a theory of capillarity. Such a theory is able to take account of systems 
in which fluid interfaces are present. The internal capillarity associated with 
fluids is the simplest case. A mathematical limit analysis associated with the 
thickness of the interface when the size of the layer goes toward zero and the 
behavior of the layer between fluid phases yield the model of material surfaces. 
Such a theory is able to calculate the superficial tension as well in the case of 
thin interfaces as thick interfaces. The static model in continuum mechanics 
of second gradient theory is extended to dynamics. The equation of motion is 
able to induce a stress tensor. In the fluid case, the theory does not lead to an 
isotropic stress tensor. 

Many developments have been done since paper 144)1 by the author: 

It is possible to obtain the radius of nucleation of microscopic drops and bub- 
bles and to develop a macroscopic theory as the Laplace theory of capillarity 
for curved interfaces [32]- The stability of interfaces is investigated with dif- 
ferential or partial derivative equations [36]. Classification of fluids endowed 
with internal capillarity is the extension of the classical fluid one's [17] . 
In fact, contact forces are of a different nature than the ones associated with 
the Cauchy stress tensor. Classical conditions with the tetrahedron construc- 
tion due to Cauchy are not efficient to study the non-linear behavior of fluids 
endowed with internal capillarity. We deduced contact forces concentrated 
on edges representing the boundaries of surfaces of separation. For example, 
such conditions are necessary to study the stability of thin films in contact 
with solid walls and the connection with the mean field molecular theories 
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The model dynamically interprets phenomena in the vicinity of the critical 
point and authorizes an investigation - at least qualitative - of the dynamic 
change of phases between bulks in fluids or fluids mixtures - New math- 

ematical equations associated with the hyperbolicity in continuous media can 
be considered in second gradient theory and consequently in internal capillar- 
ity [55|56] . 




More recently, we must notice that fluids endowed with internal capillarity are 
able to model fluid layers at a scale of some nanometers and to recognize the 
disjoining pressure concept for very thin liquid films [52] . The model can be 
applied in vegetal biology to interpret the ascent of the crude sap in very high 
trees as sequoia and giant eucalyptus [5Tf58] 
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